Pairing mechanism for high temperature superconductivity in the cuprates: what can we learn 

from the two-dimensional t - J model? 
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More than twenty years have passed since high temperature superconductivity in the copper oxides (cuprates) 
was discovered by J.G. Bednorz and K.A. Miiller in 1986 jl|]. Although intense theoretical and experimental 
elforts have been devoted to the investigation of this fascinating class of materials, the pairing mechanism re- 
sponsible for unprecedented high transition temperatures remains elusive. Theoretically, the difficulty lies 
in the fact that this class of materials, as doped Mott-Hubbard insulators involve strong electronic correla- 
tions, which renders conventional theoretical approaches unreliable. Recent progress in numerical simulations 
of strongly correlated electron systems in the context of tensor network representations |3, 4] makes it possible 
to get access to information encoded in the ground-state wave functions of the two-dimensional t - J model-a 
minimal model, as widely believed, to understand electronic properties of doped Mott-Hubbard insulators 0- 
[^. In this regard, an intriguing question is whether or not the two-dimensional t - J model holds the key to 
understanding high temperature superconductivity in the cuprates. As it turns out, such a key lies in a supercon- 
ducting state with mixed spin-singlet d + i-wave and spin-triplet p,.(p,,)-wave symmetries in the presence of 
an anti-ferromagnetic background |9]. Here, the d + .s-wave component in the spin-singlet channel breaks U{\) 
symmetry in the charge sector, whereas both the anti-ferromagnetic order and the spin-triplet pv(py)-wave com- 
ponent breaks 51/(2) symmetry in the spin sector. Therefore, four gapless Goldstone modes occur. However, 
even if we resort to the Kosterlitz-Thouless transition floll . only the tZ -I- j-wave superconducting component sur- 
vives thermal fluctuations. This turns three gapless Goldstone modes, arising from S U(l) symmetry breaking, 
into two degenerate soft modes, with twice the spin-triplet p^(py)-wave superconducting energy gap as their 
characteristic energy scale: one is a spin- triplet mode observed as a spin resonance mode in inelastic neutron 
scattering, the other is a spin-singlet mode observed as a Aig peak in electronic Raman scattering. The scenario 
allows us to predict that pairing is of -I- j-wave symmetry, with the two degenerate soft modes as the long- 
sought key ingredients in determining the transition temperature T^, thus offering a possible way to resolve the 
controversy regarding the elusive mechanism for high temperature superconductivity in the cuprates. 

PACS numbers: 74.20.-z, 74.20.Mn, 74.20.Rp 



Imagine if we would have been able to solve a model 
system describing doped Mott-Hubbard insulators on a two- 
dimensional square lattice, whose ground-state wave function 
is a superconducting state with mixed spin-singlet d + i-wave 
and spin-triplet p v(/?v)-wave symmetries in the presence of an 
anti-ferromagnetic background, with the order parameters for 
the i-wave, t/-wave, and pxip^)-wave superconducting com- 
ponents, together with the anti-ferromagnetic order parame- 
ter, shown in Fig. [1] in a proper doping range. Note that 
and As are, respectively, the spin-singlet li-wave and i-wave 
superconducting energy gaps, whereas is the spin-triplet 
;'A(/'v)-wave superconducting energy gap and is the anti- 
ferromagnetic Neel order parameter. A few peculiar features 
of this state are: (i) Both the 90 degree (four-fold) rotation 
symmetry and the translation symmetry under one-site shifts 
are spontaneously broken on the square lattice, (ii) Spin- 
rotation symmetry SU(2) is spontaneously broken, due to 
the simultaneous occurrence of both the p^{py)-wa\e super- 
conducting component and the anti-feiTomagnetic order, (iii) 
t/(l) symmetry in the charge sector is spontaneously broken, 
due to pairing in both spin-singlet and spin-triplet channels. 
Here, we emphasize that the symmetry mixing of the spin- 
singlet and spin-triplet channels arises from the spin-rotation 
symmetry breaking, simply because spin is not a good quan- 
tum number, (iv) All superconducting components are homo- 
geneous, in the sense that their superconducting order param- 



eters are independent of sites on the lattice. 

Now let us switch on thermal fluctuations. Suppose we 
restrict ourselves to a strict two dimensional system. Then, 
even if the Kosterlitz-Thouless transition jiotl is invoked, only 
spin-singlet d + s-wave superconducting component survives 
thermal fluctuations. However, the non-abelian SU{2) sym- 
metry is not allowed to be broken at any finite tempera- 
ture I n, 13]. This immediately implies that the Goldstone 
modes arising from the spontaneous symmetry breaking of 
S U(2) in the spin sector have to be turned into degenerate 
soft modes, with twice the spin-triplet px{pv)-^a^t supercon- 
ducting energy gap as their characteristic energy scale: one is 
a spin-triplet mode associated with the anti-ferromagnetic or- 
der, with the momentum transfer {n, n), and the other is a spin- 
singlet mode associated with the spin-triplet p^(py)-wave su- 
perconducting component, with the momentum transfer (0, 0). 
On the other hand, there is nothing to prevent from the break- 
ing of the discrete four-fold rotation symmetry on the square 
lattice. Actually, this broken symmetry not only manifests it- 
self in the admixture of a small i-wave component to the dom- 
inant li-wave superconducting state (see Fig.[T] left panel), but 
also protects the spin-singlet soft mode that is unidirectional 
as it arises from the pjc(py)-wave superconducting component. 

Our argument leads to a scenario that, at any finite tem- 
perature, the pairing is of li -t- .s-wave symmetry, with two 
degenerate soft modes acting as the key ingredients in deter- 
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mining the transition temperature Tc- Actually, two distinct 
energy scales 2A* and E^^s are involved, in a marked contrast 
with the conventional superconductors: 2A* arises from the 
anti-ferromagnetic Neel order parameter A^, which is respon- 
sible for pairing, with its coupling strength decreasing almost 
linearly with doping, whereas iSies - 2Ap, which is respon- 
sible for condensation. Therefore, E^^s must scale with the 
superconducting transition temperature Tc, i.e., Eres ~ ksTc, 
with ks being the Boltzmann constant [see Fig.[Tl right panel]. 
Similarly, 2A* scales as 2A* ~ ^pj*-' with T* being the so- 
called pseudogap temperature [Il3ill4ll . In addition, one may 
expect that fires < 2Arf, simply due to the fact that the pre- 
dominant li-wave superconducting component survives ther- 
mal fluctuations. Considering that both the superconducting 
gap Arf and the transition temperature Tc characterize the su- 
perconductivity, they should track each other in the entire dop- 
ing range, implying A^ ~ ksTc- In fact, for the t - J model, 
our simulation indicates that fiies ~ 1.25A(; |@|. This in turn 
allows us to estimate a universal coefficient k ^ 5.37 in the 
scaling relation: fires - « ksTc- 

Note that the two distinct energy scales in the underdoped 
regime are split off from one single energy scale in the (heav- 
ily) overdoped regime. This naturally results in a crossover 
from the Bose-Einstein condensation (EEC) regime to the 
Bardeen-Cooper-Schrieffer (BCS) regime, as conjectured in 
Ref. [15,1, which in turn is essentially equivalent to the phase 
fluctuation picture proposed by Emery and Kivelson lfl6ll . 
However, there is an important difference: the superconduc- 
tivity weakens in the heavily underdoped regime, not only be- 
cause of the loss of phase coherence, but also because of the 
decrease of the superconducting gap Aj with underdoping. 

Now a fundamental question is whether or not such a sce- 
nario is really relevant to the high Tc problem. This brings us 
to the phenomenology of the high temperature cuprate super- 
conductors. 

First, let us focus on the two distinct energy scales 2A* 
and fires. Physically, the two distinct energy scales mea- 
sure, respectively, the pairing strength and the coherence of 
the superffuid condensate. This naturally leads to two differ- 
ent phases: one is characterized by incoherent pairing, which 
may be identified with the pseudogap phase; the other is as- 
sociated with the emergence of a coherent condensate of su- 
perconducting pairs, which may be identified with the super- 
conducting phase of li -H .s-wave symmetry [see Fig.|2l. Evi- 
dence for the two distinct energy scales was reported in angle- 
resolved photoemission spectra |[l7l - l20ll . electronic Raman 
spectra I21n24il. scanning tunneling microscopy ll25ll . c-axis 
conductivit y || 26ll. Andreev reflection 02711 . magnetic penetra- 
tion depth [/8j, and other probes (for a review, see, Ref. ll29ll . 
Actually, these studies indicate that the gap near the antinodal 
region, which is identified as the pseudogap, does not scale 
with Tc in the underdoped regime, whereas the gap near the 
nodal region may be identified as the superconducting order 
parameter in the cuprates fl 17142311 . This identification offers a 
natural explanation why the two distinct energy scales in the 
underdoped regime merge into one single energy scale in the 
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FIG. 1: (color online) The doping dependence of the order pa- 
rameters for a model system describing doped Mott-Hubbard in- 
sulators on a two-dimensional square lattice, whose ground-state 
wave function is a superconducting state with mixed spin-singlet 
d + 5- wave and spin-triplet p.v(p,0-wave symmetries in the presence 
of an anti-ferromagnetic background, with the order parameters for 
the i-wave, rf-wave, and pj-(pv)-wave superconducting components, 
together with the anti-ferromagnetic order parameter A'. Left panel: 
The superconducting gaps and Aj for the spin-singlet rf-wave and 
.s-wave components as a function of doping S, with their ratio Aj/Aj 
shown in the inset. Right panel: The energy scales 2A* ~ kgT' ~ N 
and £jes = 2Ap for the anti-ferromagnetic order and the spin-triplet 
P.v(Py)-wave superconducting order as a function of doping S, with 
A' being the anti-ferromagnetic order, and kg the Boltzmann constant. 
A crossover from the Bose-Einstein condensation (BEC) regime to 
the Bardeen-Cooper-SchriefFer (BCS) regime occurs, when the two 
energy scales merge into one single energy scale in the (heavily) 
overdoped regime. Note that, in general, E^^i < 2Arf. Indeed, 
iiies ~ 1.25Aj, as predicted from the two-dimensional f-7 model. In 
addition, ii,es scales with the superconducting transition temperature 
Tc'- £ics ~ kgTc 



(heavily) overdoped regime as a consequence of the evolution 
of the Fermi arcs in the underdoped regime to a large Fermi 
surface in the (heavily) overdoped regime 1 3^ 31 1 . We em- 
phasize that the pseudogap near the antinodal region does not 
characterize a precursor to the superconducting state, in the 
sense that the pseudogap smoothly evolves into the supercon- 
ducting gap at Tc (l3. 3^ 31 1. Instead, it coexists with the 
superconducting gap of the d + s-wave symmetry in the su- 
perconducting state. More likely, a precursor pairing occurs 
in the nodal region l33l . with its onset temperature lower than 
T*, but above Tc, which may be identified with the Nernst 
regime [33, 34]. As observed, the superconducting gap near 
the nodes scales as ksTc lf29ll . This makes a strong case for 
our argument, if one takes into account the smallness of the 
i-wave superconducting gap. On the other hand, ample evi- 
dence has been accumulated, over the years, for the universal 
scaling relation fires - k kgTc, valid for both soft modes, i.e., 
the spin-triplet resonance mode in inelastic neutron scattering 
experiments 11351 - 14 ill and the spin-singlet mode observed as a 



Aig peak in electronic Raman scattering |21, 42-4^, respec 



tively. The experimentally determined k is around 6, quite 
close to our theoretical estimate. This presents a possible res- 
olution to the mysterious A\g problem [45]: the Aig mode is 
a charge collective mode as a bound state of (quasiparticle) 
singlet pairs originating from the fluctuating p v(py)-wave su- 
perconducting order. 

Second, is the pairing symmetry really of li -H i-wave na- 
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FIG. 2: (color online) A schematic phase diagram of the hole-doped 
high temperature cuprate superconductors plotted as a function of 
temperature T and hole doping S. Here, SC, PG, and AF stand for 
the superconducting, pseudogap, and anti-ferromagnetic phases, re- 
spectively. Tc, T', and Tn represent, respectively, the superconduct- 
ing transition temperature, the pseudogap temperature, and the anti- 
ferromagnetic Neel temperature. In our scenario, the PG phase is 
characterized by incoherent preformed pairs, with a pseudogap open- 
ing near the anti-nodal region, leaving the remnant gapless Fermi 
arcs near the nodes in the momentum space. In addition, the four- 
fold rotation symmetry and the translation symmetry are broken on 
the square lattice. Superconducivity occurs, when long range phase 
coherence develops at T^. 



ture? Many theoretical proposals, mainly based on the res- 
onating valence bond scenario |2], predicted a pure t/-wave 
superconductivity in the cuprates [6, 47-49]. Although a vari- 
ety of experiments have demonstrated a predominant cZ-wave 
gap fsd-^S^, strong evidence points to an admixture of an 
5-wave component to the dominant d-wave superconductiv- 
ity in muon spin rotation studies ifssi Is^ll . electronic Raman 
measurements fs^l, angle-resolved electron tunneling exper- 
iments [56], and neutron crystal-field spectroscopy experi- 
ments 15711 . In addition, a universal scaling relation of the 
superfluid density, p.s(O), at absolute zero, with the product of 
the dc conductivity crddTc), measured at Tc, and the transi- 
tion temperature Tc, indicates that a pure t/-wave supercon- 
ductivity is realized in the cuprates pSj. This scaling relation 
may be regarded as a modified form of the Uemura relation 
between the sup erfluid density ps(0) and the transition tem- 
perature Tc i59ll60ll . which works reasonably well in the un- 



derdoped regime. However, a significant deviation from the 
scaUng relation was subsequently observed ll6lll . with a salient 
feature that the deviation increases with doping. This feature 
strongly suggests that the discrepancy should be accounted for 
by removing an extra contribution from the i-wave compo- 
nent in the context of the d + i-wave pairing symmetry. This 
issue has been addressed recently li6Zl . thus supporting our 
scenario. 

Third, any theory regarding the underlying mechanism of 
high temperature superconductivity would not be complete. 



if the stripe states (for reviews, see, e.g., Ref. 163146611 ) were 
not touched upon. Surprisingly, a stripe-like state, i.e., a state 
with charge and spin density wave order, coexisting with a 
spin-triplet /7^(/:iy)-wave superconducting state, does occur as 
a ground state in the f-7 model for dopings up tO(5 ~ 0.18 |9], 
with 7/f = 0.4. Again, all symmetries, including the four- 
fold rotation and translation lattice symmetry, S U{2) spin ro- 
tation and U(l) charge symmetry, are spontaneously broken. 
A similar line of reasoning yields that, only the charge density 
wave order survives thermal fluctuations, with the concomi- 
tant occurrence of the soft modes: they arise from the S U{2) 
symmetry breaking of the spin density wave order and the 
spin-triplet px{py)-wave superconducting order, with twice 
the spin-triplet pv(/?y)-wave superconducting energy gap as 
their characteristic energy scale. Although it remains uncer- 
tain whether or not such a commensurate stripe-like state is 
an artifact of our choice of the unit cell for the tensor net- 
work representation of quantum states, an important lesson 
we have learned from our simulation is that, the t - J model 
exhibits a strong tendency towards a stripe state in the under- 
doped regime, consistent with the density matrix renormal- 
ization group |67] for a more realistic stripe pattern at doping 
1/8 and J/t - 0.35. From this we conclude that (i) static 
charge density wave order is compatible with the supercon- 
ductivity, so its possible role is deserved to be explored in the 
formation of the pseudogap; (ii) static spin density wave or- 
der is detrimental to the d+s-wave superconductivity, because 
no d + i-wave superconducting component coexists with the 
charge and spin density order in the ground state; (iii) fluc- 
tuating spin density wave order, together with the fluctuating 
spin-triplet p^{py)-wa\e superconducting order, equally well 
account for the Bose-Einstein condensation in our scenario. 
Therefore, the fluctuating stripe order is intrinsic to many, if 
not all, famiUes of the high temperature cuprate superconduc- 
tors. 

Fourth, does our prediction about the spontaneous breaking 
of the four-fold rotation symmetry and the translation sym- 
metry under one-site shifts in the pseudogap phase represent 
a physical reality? In our scenario, the PG phase is character- 
ized by incoherent preformed pairs, which occur in the anti- 
nodal regime in the momentum space, leaving the remnant 
gapless Fermi arcs in the nodal regime. In addition, the four- 
fold rotation symmetry and the translation symmetry are bro- 
ken on the square lattice. In the superconducting phase, the 
four-fold rotation symmetry breaking manifests itself in the 
admixture of an .s-wave component to the dominant d-wave 
state. In the pseudo phase, the broken symmetry protects the 
spin-singlet soft mode that is unidirectional as it arises from 
the p^(py)-wa\e superconducting fluctuations, as well as the 
fluctuating spin density order and possibly static charge den- 
sity order One may expect that measurable physical effects 
arise from the coupling of electrons with the unidirectional 
spin-singlet soft mode. Indeed, a large in-plane anisotropy of 
the Nernst effect in YBa2Cu30y was reported that sets in ex- 
actly at the pseudogap temperature T* |68]. This broken four- 
fold rotation symmetry was also detected by resistivity It69ll 
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and inelastic neutron scattering 1 7C - T^I at low doping, and by 
scanning tunneling spectroscopy 173074] at low temperature, 
without a clear connection to the pseudogap temperature T* . 

Fifth, what is the real bosonic glue for high temperature 
superconductivity in the cuprates? Many researchers have 
suggested various candidates, such as the spin resonance 
mode 1 75, [t^, the s pin fluctuation spectrum |22-81|, the 
phonon spectrum ['82', ^83'], a mode responsible for the dip ob- 
served in angle -resolved photoemission spectra |84], and the 
anti-ferromagnetic fluctuations Issilsill . as a possible bosonic 
glue. However, our scenario unveils that the spin resonance 
and the Raman A\g modes play an analogous role to a ro- 
ton in the superfluid '^He, in contrast to the conventional low- 
temperature superconductors with phonon-mediated pairing. 
This strongly suggests that, there is no bosonic glue to pair 
electrons in the cuprates, as advocated by Anderson IstIi . That 
is, pairing is an emergent phenomena; it is resulted from the 
many-body correlations in a doped Mott-Hubbard insulator: 
both the t/-wave and the extended s-wave pairings are realized 
to avoid strong on-site Coulomb repulsion. Here, we remark 
that the spin resonance and the Raman Aig modes have been 



suggested by Uemura 118811 to be roton analogues, by invoking 
soft modes in the spin and charge channels in an incommen- 
surate stripe state (see also Ref. ll89ll for an alternative expla- 
nation). In this sense, our scenario unlocks the secret of the 
yet-to-be-achieved solid, as an analogue of the solid "^He, in 
the cuprates. 

In closing, let us discuss a possible candidate of the model 
system that is able to reproduce all the features needed to ad- 
dress the high Tc problem, as mentioned at the beginning of 
the text. With the fact in mind that the two-dimensional t - J 
model does exhibit a + i-wave superconducting state and a 
stripe-like state in a certain doing range for a physically re- 
alistic parameter regime 1^, we speculate that, most likely, 
an extension of the t - J model to include the next-nearest- 
neighbor and the third-nearest-neighbor hopping terms, which 
are known to be significant in the cuprates, will do the job. 
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